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ABSTRACT 

The origin of rovibrational H2 emission in the central galaxies of cooling flow clusters 
is poorly understood. Here we address this issue using data from our near-infrared 
spectroscopic survey of 32 of the most line-luminous such systems, presented in the 
companion paper by Edge et al. (2002). 

We consider excitation by X-rays from the surrounding intracluster medium 
(ICM), UV radiation from young stars, and shocks. The v=l-0 K-band lines with 
upper levels within 10^ K of the ground state appear to be mostly thermalised (im- 
plying gas densities > 10^ cm~^ ), with the excitation temperature typically exceeding 
2000 K, as found earlier by Jaffe, Bremer & van der Werf (2001). Together with the 
lack of strong v=2-0 lines in the H-band, this rules out UV radiative fluorescence. 

Using the CLOUDY photoionisation code, we deduce that the H2 lines can origi- 
nate in a population of dense clouds, exposed to the same hot (T ^ 50000 K) stellar 
continuum as the lower density gas which produces the bulk of the forbidden optical 
line emission in the Ha-luminous systems. This dense gas may be in the form of self- 
gravitating clouds deposited directly by the cooling flow, or may instead be produced 
in the high-pressure zones behind strong shocks. Furthermore, the shocked gas is likely 
to be gravitationally unstable, so collisions between the larger clouds may lead to the 
formation of globular clusters. 

Key virords: galaxies: active - galaxies:starburst - galaxies: cooling flow - X- 
rays:cooling flow 



1 INTRODUCTION 

In the companion paper by Edge et al. (2002) (hereafter pa- 
per 1) we presented H and K-band spectra of 32 of the most 
Hq line-luminous cooling flow (CF) central cluster galaxies 
(CCGs) in the i? 05^ T Brightest Cluster Sample (Crawford 
et al. 1999). One or more of the rovibrational H2 v=l-0 
S(0)-S(7) lines were detected in 23 systems; one or both 
of [Fell] AAl. 258, 1.644 were seen in 14 systems, and there 
is also evidence for higher excitation H2 lines (e.g. v=2-l 
S(l), S(3)) and some coronal lines ([Si VI], [Si XI], [S XI] 
and [Ca VIII]). This sample quadruples the number of CF 
CCGs with H2 detections, and builds on work by Jaffe & 
Bremer (1997), Falcke et al. (1998), Donahue et al. (2000) 
and Jaffe, Bremer & van der Werf (2001). Now that warm 
(~ 1000 — 2500 K) molecular hydrogen is known to accom- 
pany ionized material in the core of a flow, our aim in this 



paper is to shed light on the associated excitation mecha- 
nisms. 

Jaffe & Bremer (1997) found that the H2 emission in CF 
CCGs is too luminous to be simply due to material passing 
through ~ 2000 K- the temperature at which it is collision- 
ally excited - whilst cooling from ~ 10^ K (the well-known 
'H rec' problem). Based on a high v=l-0 S(l)/Ha ratio, 
they implicated coUisional excitation by suprathermal elec- 
trons liberated by X-ray photoionization deep within the 
cold, molecular core of the line-emitting clouds; the X-rays 
were hypothesized to originate in the surrounding intra- 
cluster medium (ICM). A similar model was advanced by 
Wilman et al. (2000) for Cygnus A, which exhibits H2, H 
recombination and [Fell] emission out to radii of 5kpc. Flux 
ratios involving the v=l-0 S(l), S(3) and S(5) lines deviate 
markedly from the predictions of LTE, suggesting that non- 
thermal excitation is important. If the latter is provided by 
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hard X-rays from the obscured quasar nucleus, the imphed 
H2 mass within 5 kpc of the nucleus is ~ 10^° M0 , account- 
ing for 10® yr of mass deposition from the CF. 

From their HST NICMOS narrow-band imaging of the 
K-band H2 emission in three CF CCGs, Donahue et al. de- 
duced that the emission in NGC 1275 is principally due to 
the AGN (as did Krabbe at al. 2000), but that the more 
spatially-extended emission in A2597 and PKS 0745-191 is 
probably due to UV radiation from young stars. They based 
this conclusion on the close morphological correspondence 
between the H2 and optical line emission, and on a compar- 
ison of the H2/Ha ratios with predictions for shock models. 
X-ray heating by the ICM, and UV fluorescence by stars and 
AGN. 

Most recently, Jaffe, Bremer & van der Werf (2001) pre- 
sented UKIRT CGS4 K-band spectra of 7 well-known cool- 
ing flows. They observed that the relative strengths of the 
v=l-0 lines agree well with the predictions of LTE, but that 
there are deviations amongst the higher vibrational states. 
The thermal nature of the line production implies that the 
molecular gas is overpressurised by 2 to 3 orders of mag- 
nitude with respect to the X-ray and optical line emitting 
components. They sought to explain this using an ablative 
or 'rocket' model of the gas system in which X- or UV- 
radiation heats the surface layers of a cold, gravitationally 
bound molecular cloud to 2000 K, at which temperature the 
material ceases to be gravitationally bound and expands into 
space, ionizes, and emits the optical lines. They deduced ion- 
ized to molecular line ratios which are lower than those of 
starburst regions, indicating that alternative heating mech- 
anisms are necessary. 

Of relevance to the analysis presented in this paper are 
the results of the CO survey for cold molecular gas in CF 
CCGs presented by Edge (2001), many of whose targets are 
also in our UKIRT survey. He detected CO emission in 16 
of these CFs, consistent with jvIq of molecular gas 

at ~ 40K (assuming a standard C0:H2 conversion factor). 
The implied mass of cool H2 correlates better with the Ha 
luminosity than with the global X-ray mass deposition rate 
of the CF, suggesting that young stars are warming a pop- 
ulation of molecular clouds. Since the physical conditions 
in the material probed by the K-band H2 lines are likely 
to occupy an intermediate regime between the optical and 
CO line-emitting components, they offer the potential for 
further investigation into the link between the CO and Ha 
emission. Indeed, we showed in paper 1 that there is a good 
correlation between the H2 mass inferred from the CO emis- 
sion and the 1-0 S(l) line luminosity, with the exceptions of 
the strong radio sources Cygnus A and PKS 0745-191. 



2 ANALYSIS 

We address the H2 excitation mechanism from several an- 
gles: firstly, by considering the relative strengths of the vari- 
ous H2 lines; secondly, by comparing the fluxes in the H2 and 
H recombination lines and using established results for the 
origin of the optical line emission in these systems. Through- 
out we draw upon our own and published calculations for the 
H2 emission under various physical conditions (viz. irradia- 
tion by UV and X-rays, and in shocks), as well as empirical 
results for the H2 emission in Seyfert and starburst galax- 



ies. A comparison between the widths of the H2 and CO 
lines (the latter are determined to an accuracy of 50kms~^ , 
and generally lie in the range 100-300 kms"^ FWHM; Edge 
2001) would also elucidate the link between the CO and H2 
emission, but is not possible due to the low resolution of our 
UKIRT spectra (~ 600 km s"^ FWHM at the centre of the 
K-band). 

2.1 H2 line ratios 

We begin by comparing the observed H2 line ratios with the 
predictions for thermal (coUisional) excitation, which domi- 
nates in molecular gas with a density riT ^ 10* cm"'' , heated 
to a few thousand Kelvin by a source of X or UV radia- 
tion, or in shocks. Under these conditions, the occupation 
numbers of the excited rotation-vibration levels of the H2 
molecule will be in thermal equilibrium at a temperature Tex 
equal to the kinetic temperature of the gas. Therefore, for a 
given object, the flux in the ith emission line will satisfy 
the relation logi^iXi / AiQi) = constant — Ti/Tex, where Ai is 
the wavelength of the line, and Ai, Ti and are the sponta- 
neous emission coefficient, excitation temperature and sta- 
tistical weight of the upper level of the transition (the latter 
assumes an ortho:para H2 abundance ratio of 3:1, appropri- 
ate where collisions dominate). LTE plots were constructed 
for 14 objects (those in which 3 or more H2 lines were de- 
tected), and are shown in Figs. |l| and ^ for the September 
1999 and March 2000 objects, respectively. 

These figures show that where reliable measurement is 
possible, the v=l-0 lines with upper levels within ~ lO** K of 
the ground state are thermally excited at temperatures in 
excess of 1600 K, with most being above 2000 K. The only 
exception to this is All, where the excitation temperature 
for the v=l-0 S(l), S(2) and S(3) lines is 890 ± 200 K. Im- 
plied temperatures above 4000 K are unphysical since H2 
is rapidly dissociated by collisions in hotter gas. Thus the 
result for A1795 should be regarded as spurious, especially 
since it conflicts with the results of Jaffe et al. (2001) who 
found line ratios in LTE at approximately 2200 K. The dis- 
crepancy is probably due to our spectrum having a lower 
signal to noise ratio and the difficulty of assessing the true 
continuum underlying the v=l-0 S(3) and S(5) lines which 
fall in a region of low atmospheric transmission. The ther- 
malised nature of the lower-lying transitions thus demands 
that the total atomic and molecular hydrogen density exceed 
~ 10* cm~^ , the critical density for coUisional deexcitation. 
In combination with the measured temperatures, the implied 
pressures exceed those in the X-ray and optical line emit- 
ting systems by 2 to 3 orders of magnitude, and we therefore 
confirm the result of Jaffe et al. (2001). Together with the 
absence of prominent v=2-0 lines in our H-band spectra, 
this also rules out UV radiative fluorescence in cold, low 
density gas. The latter mechanism was advocated by Don- 
ahue et al. for A2597 and PKS 0745-191 based on the value 
of the v=l-0 S(3)/Ha ratio and energetic arguments. UV 
fluorescence may, however, play some role in the production 
of the higher excitation lines such as v=l-0 S(7) and v=2-l 
S(l) and S(3) which are seen to deviate from LTE in Figs. 1 
and 2. CoUisional excitation by a suprathermal photoelec- 
tron (as modelled by Maloney, HoUenbach & Tielens 1996) 
is an alternative non-thermal excitation mechanism which 
may also supply some of the flux in these lines. 
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Figure 1. LTE diagnostic plots for objects in the September 1999 dataset with three or more H2 line detections. 
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Once in the thermal regime, it is difficult to deter- 
mine the heat source for the gas from the H2 lines alone. 
Mouri (1994) has, however, shown that there is a demar- 
cation in excitation temperature between heating by UV 
on the one hand, and by shocks or X-rays on the other; 
from the models of Sternberg & Dalgarno (1989) (hereafter 
SD89), he showed that UV radiation can heat the molecular 
gas to no higher than 1000 K, whilst X-ray heating (from 
the models of Lepp & McCray 1983) produces tempera- 
tures ~ 2000 K; from the models of Brand et al. (1989), 
shocks occupy the upper end of the range in between. Thus, 
the temperatures derived from Figs.^ and |2| appear to rule 
out UV heating and marginally favour X-ray heating over 
shocks. There are two caveats to this conclusion. Firstly, X- 
ray heating has been shown to be energetically incapable 
of powering the optical line emission in CF CCGs (see e.g. 
Crawford & Fabian 1992). Secondly, many of the most lu- 
minous systems (with L(Hq) > lO^^ergs"^, i.e. most of 
the objects in our sample) have substantial populations of 
young stars, which produce enough UV radiation to power 
the observed Ha emission (Crawford et al. 1999). The op- 
tical emission line spectra differ from those of 'normal' Ifll 
regions, with strong Balmer lines relative to forbidden lines 
of [OUT], [Nil] and [SIX]; they can be reproduced with very 
hot stars {T > 40000 K) and a low mean ionisation param- 
eter U in the range -3.5 to -4.0 (Allen 1995). 

2.2 H2 and HII line emission 

The HST NICMOS emission line imaging study of Donahue 
et al. revealed a close spatial correspondence between the H2 
and optical Ho? emission, suggesting a common excitation 
mechanism. 

We begin with an empirical comparison of the H2 and 
HII line emission. Jaffe et al. (2001) noted that in Seyfert, 
starburst and ultraluminous infrared galaxies (ULIRGs), the 
ratio v=l-0 S(l)/Br7 < 1. Assuming case B recombination 
this implies that v=l-0 S(l)/PaQf < 0.083 for such objects. 
Fig. ^ shows, however, that for all our targets with de- 
tectable S(l), v=l-0 S(l)/PaQ > 0.2. Thus, whatever the re- 
lationship between the H2 and HII line emission in Seyferts, 
starburst galaxies and ULIRGs - and this is by no means 
certain in itself - it does not apply to our CF galaxies. Fig. ^ 
also shows no correlation with the power of the central ra- 
dio source, which might have been expected if the associated 
AGN preferentially enhanced the production of the HII or 
H2 lines. 

Following Donahue et al. we also compare the observed 
v=l-0 S(l)/Ha ratio with the predictions of various models, 
as shown in Fig. ^ Since the Ha and S(l) fluxes are both 
slit fluxes derived from separate observations at different 
position angles, this may introduce a scatter as large as a 
factor of 4 in v=l-0 S(1)/Hq. This is because the emission 
line images in Donahue et al. show that as much as half of 
the H2 emission (and more of the Ha) arises in filamentary 
structures at radii greater than 1 arcsecond. 

For UV excitation the v=l-0 S(l)/Ha ratio is sensi- 
tive to the shape of the exciting spectrum near the Lyman 
edge: photons with A < 912A produce Ha, whilst those with 
912 < A < llOSA (i.e. in the Lyman and Werner bands) 
lead to H2 line emission. For a power-law spectrum typi- 
cal of an AGN, they infer from the results of Black & van 
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Figure 3. v=l-0 S(l)/Paa versus monochromatic 1.4 GHz 
radio power. 

Dishoeck (1987) that the S(l)/Ha ratio lies in the range 
0.0033-0.017. Introducing a break by a factor 100 at the 
Lyman edge, to mimic a UV radiation field dominated by 
stars, yields a higher ratio of S(l)/Ha ~ 0.33 — 2. It should 
be noted, however, that the Black & van Dishoeck models 
do not include coUisional excitation processes for the vibra- 
tionally excited states and are thus not applicable to densi- 
ties above 10^ cm~^ . 

Donahue et al. also used the models of Maloney, Hol- 
lenbach & Tielens (1996) to compute the S(l)/Ha ratio 
expected for X-ray heating by the ICM. For a thermal 
bremsstralung spectrum with T = 8 keV and a constant 
gas pressure of 10* cm""^ K, S(l)/Ha varies from 0.0125 to 
0.125 as the total column density of the cloud is increased 
from 5 X 10^^ to 1.5 x 10^^ cm . Since hard X-rays can 
excite H2 emission through collision with suprathermal pho- 
toelectrons, the S(l)/Ha can be increased further by aug- 
menting the column density. Uncertainties in the shape of 
the UV/soft X-ray spectrum can also affect the contribution 
from the highly ionised layer at the cloud surface. Not with- 
standing such uncertainties. Fig. ^ shows that of the mech- 
anisms considered by Donahue et al. X-ray heating is the 
most viable for the bulk of the objects, as found in section 
2.1. However, in the next subsection we present calculations 
for thermal excitation in dense gas heated by the UV radia- 
tion from young stars, and describe how it provides a more 
natural explanation for the H2 emission in these systems. 

2.3 CLOUDY photoionisation modelling 

To model thermal excitation of the H2 lines in dense gas 
heated by UV radiation (as opposed to the non-thermal 
UV fluorescence modelled by Black & van Dishoeck), we 
use the photoionisation code CLOUDY (version 94.00; Fer- 
land 1996). As mentioned in section 2.1, we adopt a Kurucz 
model atmosphere at 50000 K, with the ionisation parame- 
ter U in the range -3.5 to -4.0 when incident on gas with 
hydrogen number density ri{H) = 200 cm"'' , in order to 
reproduce the optical emission line spectrum (this density 
matches that indicated by the [SII]AA6717, 6731 doublet, see 
e.g. Voit & Donahue 1997 for A2597). Note that CLOUDY 
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Figure 4. v=l-0 S(l)/Ha versus Ha luminosity where 
known (mostly taken from Crawford et al. 1999, corrected 
for intrinsic reddening which is typically E{B — V) ^ 0.3). 
The regions bounded by the pairs of solid and dotted lines 

correspond, according to Donahue ct al., to UV fluorescence 
by an AGN and stellar sources, respectively. The dashed 
lines encompass the range for X-ray heating by the ICM. 
As discussed in section 2.3, however, thermal excitation by 
stellar UV radiation provides the most natural explanation 
for both the Ha and H2 emission. 



does not list the H2 lines separately but instead gives the 
'sum of the fluxes of the H2 lines near 2/xm' (which we de- 
note by 2/^2). Since v=l-0 S(l) is typically the strongest H2 
line in our spectra, we do not for the moment draw a distinc- 
tion between Ei/2 and v=l-0 S(l) in comparing model with 
observation. In the models we neglect dust and any contri- 
bution from magnetic pressure. Later in this sub-section we 
briefly discuss the effects of dust on the models and estimate 
the actual value of the 1-0 S(l)/E//2 ratio. 

The inferred pressure difference between the optical and 
infrared emission line components clearly means that they 
cannot come from different parts of a single cloud in pressure 
equilibrium with the X-ray cluster gas. Such clouds would 
in any case have T,H2/Ha ratios several orders of magni- 
tude below the observed range shown in Fig. 4. We can also 
dismiss the possibility that the two components arise at dif- 
ferent depths within the same self-gravitating clouds, with 
the optical lines being produced in the surface layers which 
are in pressure balance with the diffuse X-ray gas: assuming 
hydrostatic equilibrium, the masses of the individual clouds 
would have to exceed 5 x W^Rp^ Mq in order to produce the 
requisite pressure difference between the two components 
{Rpc being the cloud radius in pc). 

A more likely possibility is that the H2 lines originate 
in a separate population of much denser, self-gravitating, 
clouds. As an example, consider clouds with n{H) = 
10'' cm^'^ exposed to the above stellar continuum with 
U = —6.2 (i.e. with the same incident flux as that seen 
by low density, constant pressure clouds with n{H) = 
200 cm~* and U = —3.5). The resulting emission line 
spectrum has E_ff2/Ha — 0.1 and the optical forbid- 
den lines from ionised species are coUisionally suppressed 
viz. [NII]A6584/Ha = 0.042, [SII]A6717/Ha = 0.033, but 
[OI]A6300/Ha = 0.12. Furthermore, the kinetic temperar 



ture in the H2-emitting gas is 1000-2000 K and the pressure 
is several times 10* cm~^ K, as required. 

We thus propose that the variations in the [Nil] /Ha, 
[SII]/Ha and [OIII]/H/3 ratios with Ha luminosity, as shown 
in Figs. 7 and 8 of Crawford ct al. (1999), and likewise the 
variation of Ei/2/Ha shown in Fig. 4 above, can be repro- 
duced by a mixture of three separate emission line compo- 
nents. The relative line intensities within each component 
are fixed, but the fraction of the total Ha emission that 
each contributes varies with the Ha luminosity. The three 
components are: 

(a) Some mechanism (e.g. mixing layers; Crawford & 
Fabian 1992) which produces the emission line ratios char- 
acteristic of the 'type I' systems with Ha luminosities be- 
low 10"^ erg s"^; A2199 and A2052 arc good examples of 
such systems, having [NII]/Ha ~ 2.3, [SII]/Ha ~ 0.7, 
[0III]/H/3 ~ 1 (but with a large scatter), and no appre- 
ciable H2 emission. 

(b) The spectrum produced by low density gas {n{H) ~ 
200 cm~^ ) in pressure equilibrium with the X-ray cluster 
gas, when photoionised by the continuum of hot massive 
stars with U in the range -3.5 to -4. As discussed, this re- 
produces the emission line ratios of the luminous 'type H' 
systems, with [NII]/Ha ~ 0.6, [SII]/Ha ~ 0.3, [0III]/H/3 ~ 
0.75, and again no appreciable H2 emission. 

(c) A population of dense {n{H) > 10* - lO^cm'^), 
perhaps self-gravitating, clouds which are exposed to the 
same photoionising flux as those in (b); these produce the 
H2 lines, with Ef/'2/Ha ~ 0.03 - 0.50, and very little flux in 
the forbidden hues of [Nil], [SII] and [OIII]. 

Suppose that the fractions of the Ha luminosity con- 
tributed by components (a), (b) and (c) are /a, fb and /c, 
respectively, and that /b = /c = 0.5(1 — /a). Now if /a 
decreases with Ha luminosity from essentially unity at log 
L(Ha)=40 to around 0.2 at log L(Ha)=42, with a functional 
dependence chosen to yield the variation of [Nil] /Ha with 
log L(Ha) in Fig. 8 of Crawford et al. (1999), then this 
also reproduces the variations of the [SII]/Ha, [OIII]/H/3 
and E/i'2/Ha ratios with luminosity. In particular, the ab- 
sence of H2 emission in the lower luminosity objects (log 
L(Ha)< 41) is accounted for. This phenomenological result, 
in which two or more components of emission appear to be 
relatively evenly-balanced over a range of more than two or- 
ders of magnitude in Ha luminosity, perhaps implies that 
some sort of feedback is operating. However, the model pa- 
rameter space and the scatter in the observed correlations 
arc sufficiently large that it is difficult to say anything more 
quantitative. The most important result to emerge from the 
CLOUDY modelling is that a separate population of dense 
clouds is required to produce the H2 emission. In the next 
subsection we discuss the origin of these clouds. 

Including dust in the CLOUDY models by means of 
the default 'grains' option increases the total H2 emission in 
the dense (component c) clouds by 20-30 per cent; the in- 
trinsic optical emission produced by the clouds is affected 
much less. Although the extinction from the illuminated 
cloud surface to the point where the model integration stops 
(the point where the temperature has fallen to 200 K), is 
negligible {Ay ~ 10~^mag), predicting the observed spec- 
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trum is not straightforward. This is because grains at greater 
depth within the molecular cloud can scatter photons back 
with a wavelength-dependent albedo: the result is highly 
geometry-dependent. It should be emphasised that the gas- 
phase abundances wore left unchanged at their solar values 
whilst the grains were introduced, so this illustrates that 
the grains per se do not have much effect. Using instead the 
'abundances ism' option in CLOUDY (which assumes deple- 
tion onto grains in accordance with observations), results in 
a reduction in the intrinsic H2/Ha ratio by a factor of ^ 6, 
owing to the existence of a much less extensive zone of warm 
(> lOOOK) H2 in this case. The elemental abundances in the 
cold clouds in cooling flows may, however, be quite different 
from those in the ISM; emission line modelling (see Allen 
1995 and references therein) suggests metallicities of around 
solar; furthermore, recent XMM-Newton grating spectra of 
several clusters show a lack of emission lines from gas cool- 
ing below 1-2 keV, one explanation for which is an inhomo- 
geneous distribution of metals within the ICM (Fabian et 
al. 2001). 

To estimate the v=l-0 S(l) flux produced by the model 
for comparison with CLOUDY's T,H2, we first extract the HI 
and H2 densities and the temperature as functions of depth. 
For the gas densities of interest, we can assume coUisional 
excitation of the H2 and use the approximations described 
by HoUenbach & Shull (1979) for computing the fractions 
of H2 in the various excited levels. Firstly, the fractions of 
molecules in each of v=0, 1 and 2 are estimated by approxi- 
mating the molecule as a three level system. Assuming LTE 
for the rotational levels within v=l, and an orthorpara ratio 
of 3:1, the fraction of molecules in the v=l, J=3 state and 
hence the v=l-0 S(l) flux follow. For the above models we 
find v=l-0 S(l)/S//2 = 0.22, close to the figure estimated 
from the observed line fiuxes listed in paper 1. To reproduce 
the 1-0 S(l)/Ha ratios in Fig. 4, we thus need component (c) 
clouds with Sff2/Ha ratios (and hence densities), towards 
the upper end of the range considered above. 

Using the same formalism, we can also calculate the 
populations of the v=0 rotational levels (which follow LTE 
up to some value J = Jm&yi, and then decrease more rapidly 
at higher J as the radiative decay rates increase), and 
hence the fluxes of the pure rotational 0-0 S series lines. 
For the component (c) clouds above, we find: 0-0 S(3)/l- 
S(l)= 0.4 ± 0.1, 0-0 S(2)/l-0 S(l)~ 0.06 and 0-0 S(l)/1- 
S(l)~ 0.04 (for our targets these are the best-placed lines, 
lying in the mid- infrared N and Q-bands). Although these 
0-0 S line fluxes are around 10 times higher than predicted 
from an LTE extrapolation of the 1-0 S series line fluxes 
(using their observed excitation temperatures of ~ 2000 K), 
they will still challenge the capabilities of forthcoming mid- 
infrared spectrographs, such as Michelle on the UKIRT and 
Gemini telescopes. A measurement of the 0-0 S(3)/l-0 S(l) 
ratio would, however, provide a useful constraint on the vari- 
ation of temperature with column density over the region of 
the cloud where the transition from atomic to molecular hy- 
drogen occurs. 

2.4 The physical picture 

One interpretation is that these dense gas clouds cooled di- 
rectly from X-ray-emitting temperatures, as posited in the 
cooling flow scenario (see Ferland, Fabian & Johnstone 1994 



for a detailed discussion of the physical conditions within 
such clouds, albeit those at larger radii in the flow). If 
the high pressures in the H2 gas are indeed due to self- 
gravitation, this is no surprise: self-gravitating molecular 
cloud cores must bo abundant in copiously star-forming en- 
vironments such as these. In this picture, the gas which pro- 
duces the forbidden optical lines (component b above), is 
that which 'boils off' these self-gravitating clouds, as pro- 
posed by Jaffe et al. (2001). 

The total mass of warm H2 in these clouds, assuming 
it to be in LTE at 2000 K, is 4.2 x 10''l/(S'l)4o M© where 
L(<S'l)4o is the luminosity of the v=l-0 S(l) line in units 
of 10"° erg s"^ (following Scoville et al. 1982). For the lu- 
minosities of interest this equates to masses of 10^^" M©. 
The gravitational collapse timescale of the clouds is tgrav = 

where ns is the particle number 
density in units of 10® cm~* . Thus, the warm molecular gas 
needs to bo roplonised at the rate of a few solar masses per 
year, a rate which is within an order of magnitude of the 
likely X-ray mass cooling rate in the central regions and the 
star-formation rate therein. 

As an alternative to self-gravity for producing the dense, 
high pressure clouds, one could appeal to shocks. Whilst 
there are considerable difficulties - discussed in detail by 
Jaffe et al. (2001) - in constructing shock models which 
can themselves (i.e. through the shock emission) produce 
the observed H2/Ha ratios and the pressure differences in a 
manner consistent with the observed kinematics, one could 
instead invoke shocks just to create the dense, high pres- 
sure clouds. This could happen as follows. Suppose that two 
of the low density clouds of component (b) above, collide 
at a relative speed of a few hundred kms^^ (i.e. equal to 
the observed velocity dispersion in the CCG). The shocks 
are strong so the density jumps by a factor of 4 through 
the shock to ~ 10^ cm~^ , and the immediate post-shock 
temperature is 10® — lO'^ K; the pressure in the post-shock 
region is then 10* — lO^'' cm""^ K, comparable to that of the 
H2-emitting region. The gas then radiatively cools and its 
density simultaneously increases to maintain the high pres- 
sure. If it remains in the isobaric post-shock region until its 
temperature has dropped to ~ 2000 K, it will then come into 
equilibrium with the ambient radiation field of the young 
stellar population and emit the observed H2 spectrum, ex- 
actly as in the case where the clouds are self-gravitating. 
The viability of this model depends on two conditions being 
satisfied: (i) on the shock geometry being such that the gas 
remains within the high pressure region until it has cooled 
to ~ 2000 K; (ii) on the line emission produced by the shocks 
themselves being negligible in comparison with that due to 
photoionisation of both the unshocked low density clouds 
and the shocked dense clouds. 

It is not clear whether either condition can be met. 
Firstly, if the clouds are magnetised, this will limit the 
post:preshock compression ratio to a maximum of rim/wo = 
77wioo/6, where wioolOOkms"^ is the shock velocity and 
the magnetic field is assumed to scale with the preshock 
density as Bo = bn%'^ jiG and 6 is of order unity (see HoUen- 
bach & McKee 1979 for the derivation of this result). Once 
the gas has reached this density, any further cooling will 
be isochoric, not isobaric, so from a pre-shock density no = 
200 cm~^ the required densities of 10® cm~^ or more will not 
be attained unless the shock speed is at least 650kms~^ . 
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Secondly, at the speeds of interest the postshock material 
will itself be a powerful emitter of UV radiation and optical 
line emission: according to Dopita & Sutherland (1995), the 
UV and Balmer line fluxes emitted by the shock front are 
8.9x 10"^Un)o'*no and 7.4 x IO^^'OiooVk) erg cm"^ s"-^ , respec- 
tively. These figures are comparable to the external stellar 
flux incident on, and the resulting H/3 line flux from, the 
dense post-shock gas in the CLOUDY model. There may 
thus be no need to invoke the existence of an additional 
stellar component: the radiation emitted by the shocked gas 
(which, for these shock speeds has a temperature of sev- 
eral times 10'' K- comparable to the temperatures of the 
massive stars invoked above) could serve the same purpose. 
This process, in which the UV radiation field and the dense 
gas are created together, could constitute the feedback re- 
ferred to in section 2.3. Concerning the H2 emission from 
the shocks themselves, we refer to the work by HoUenbach & 
McKee (1989) who simulated fast shocks (30-150 kms~^ ), in 
gas at densities no = lO"* — 10® cm~^ : such shocks are likely 
to be dissociative J-shocks from which the H2 emission is due 
to formation pumping with a spectrum which resembles that 
of UV fluorescence. For no ^ 10® cm~^ , there is a contribu- 
tion from coUisionally excited H2 in gas at ~ 500 K, but this 
contribution is only significant for the pure rotational lines 
(0-0)8(0-5). The v=l-0 S(l)/Ha ratios from such shocks lie 
below the observed range in Fig. 4 for all but the slowest 
shocks {v < 50kms^^ ). 

Shocks have previously been invoked to account for 
the optical line emission in these systems. Crawford & 
Fabian (1992) used mixing layers to account for the spec- 
tra of the 'type 1' systems (component a above), combined 
with an increasing contribution of line emission from shocks 
in dense gas (no > 10® cm~^ ) to explain the transition to 
the more Ha-luminous 'type IP systems. The latter authors 
also found correlations between the velocity dispersion and 
the [Nil] /Ha ratio, and between the local velocity disper- 
sion and the Ha surface brightness. Both correlations are 
naturally accounted for by the shock model, as discussed by 
Dopita & Sutherland (1995). 

Further insight into the shock model, and an estimate 
of the ratio of the mass emitting Ha to H2, may be ob- 
tained as follows. Suppose that within a region of radius 
R there exists a mean number density A/" of clouds of ra- 
dius r, velocity v and coUisional cross-section a ~ 7rr^. The 
rate of collisions per unit volume is Af^av, so the mass of 
shocked gas at high pressure, given that the high pressure re- 
gion lasts for a time t ~ r/v, is Mp ~ ^■KF^nm-p^irr^M^crr 
(where n is the particle density in the unshocked clouds, 
which we take to be ~ 200 crri^'^ , as found for the compo- 
nent (b) clouds, and rrip is the proton mass). The total mass 
of clouds is Mc ~ ^-KR^Nnmp^Tir^. Thus Mp/Mc = A/'vrr^ 
We may replace Mp by the mass of warm H2 which is 
given above as 4.2 x 10*I/(5'l)4o M©; similarly, we may set 
Mc equal to the mass in the Hq clouds, which is given by 
1.2 X 10®L(//a)4o Mq (using the expression for the Ha emis- 
sivity for case B recombination at 10* K in Osterbrock 1989). 
Thus Mp/Mc = 0.35L(S'l)/L(Ha) = 0.0035 - 0.035, since 
Fig. 4 shows that L(Sl)/L(Ha) generally lies in the range 
0.01-0.1 for L{Ha) > lO"*^ ergs"^ . Thus the filling factor of 
the clouds, / = ^Trr^Af, lies in the range 0.005-0.05. The 
fact that the amount of shocked gas oc JV^ might explain 



why the H2 emission is peaked towards the centre of the 
galaxy where A/" is likely to be highest. 

It is also worth noting that the shocked gas may be- 
come self-gravitating. The post-shock cooling time tcooi < 
lO'* yr (with equality holding for 10^ K), so tcooi Across 
unless r < 1 pc. The Jeans mass of the postshock gas is 
Mj = Jlfl^^.i ~ 6 X lO«T4V:24Ps'''^M0, where At-24 is 
the mean mass per particle in units of lO"^** g. It is thus 
possible that collisions between the larger clouds (r > 10 pc, 
say) can produce self-gravitating objects with masses com- 
parable to those of globular clusters, and may thus be rel- 
evant to the formation of the latter. Indeed, the specific 
frequency of globulars is sometimes high in cD galaxies (see 
e.g. the review of Harris 1991) and a young population has 
been seen at the centre of NGC 1275 (Holtzman et al. 1992). 

Since the shocked gas can become unstable to gravita- 
tional collapse and hence form stars, the gas cycling pro- 
cesses are undoubtedly complex. The above distinction be- 
tween the shock and self-gravity scenarios for producing the 
dense gas may thus be too simplistic: some of the low density 
clouds may form by evaporation from the self-gravitating 
clouds, as well as by direct cooling out of the hot ICM, and 
both stars and UV radiation from shocked gas may con- 
tribute to the excitation. 



3 CONCLUSIONS AND FUTURE PROSPECTS 

Our UKIRT spectra demonstrate that the central galax- 
ies of cooling flow clusters with Ha luminosities above 
lO'^^ergs"^ exhibit rovibrational H2 line emission at the 
level of 0.01-0.1 Ha. The relative strengths of the H2 lines 
imply that the emission is thermally excited in dense gas 
(n > 10® cm~^ ) at temperatures ~ 2000 K, and is thus over- 
pressurised by 2-3 orders of magnitude with respect to the 
optical emission line gas and the X-ray ICM. The emission 
could originate in a population of dense clouds heated by 
the young stellar populations which are known to exist in 
such systems and which can reproduce the optical forbid- 
den line spectrum when incident on lower density gas. These 
dense clouds may be self-gravitating or confined in the high- 
pressure regions behind strong shocks. The fact that the high 
and low density clouds ultimately have a similar origin and 
are excited by the same radiation field accounts for the sim- 
ilar morphologies of the H2 and optical emission line images 
found by Donahue et al. 

Future observations at higher spectral resolution could 
address the kinematics of the H2 emission for comparison 
with those of the CO emission (most of whose line widths 
arc known to within 50kms~^ and lie in the range 100- 
300 km s"'^ ; Edge 2001), which may originate at a greater 
depth within the same dense clouds. Deeper spectra of an en- 
larged sample of low Ha luminosity systems (< lO'*^ erg s~^ ) 
should also be obtained to establish at what level they emit 
H2. Fig. 4 tentatively suggests that any H2 in these sys- 
tems may be produced by a different mechanism. Our mod- 
els also open the possibility of detecting some of the purely 
rotational H2 lines in the mid-infrared with forthcoming in- 
struments such as Michelle on the UKIRT and Gemini tele- 
scopes. Observations of the 0-0 S(3)/l-0 S(l) ratio, for ex- 
ample, would provide a useful constraint on the variation 
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of temperature with depth over the HI-H2 transition region 
within the clouds. 
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